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1.  INTRODUCTION 

The  purpose  of  this  research  was  to  investigate  some 
of  the  effects  of  boundary  wall  perturbations  on  finite 
amplitude  standing  waves  in  a  rigid-walled  rectangular 
cavity. The  investigation  was  prompted  by  an  examination 
of  the  experimental  results  of  Coppens  and  Sanders [3^, 
the  research  of  DeVall[5j  and  of  Kilmer[M-]  , which  suggested 
the  existence  of  the  excitation  of  modes  other  than 
those  belonging  to  the  family  of  the  driven  mode. 


2. BACKGROUND 

A  plane  elastic  wave  travelling  in  a  non-dissipative 
fluid  will  change  waveform  as  predicted  by  the  relevant 
hydrodynamic  equations  [6j,[7j.If  the  problem  is  extended 
to  absorptive  media, only  waves  of  relatively  high  amplitude 
will  change  waveforms  appreciably. 

At  the  Naval  Postgraduate  School ,Coppens  and  Sanders 
[3], Kilmer  [t] ,and  DeVall  [5]  have  dealt  with  the  study 
of  finite  amplitude  waves  in  rigid-walled  rectangular  cavities 

One  interesting  result  of  these  cavity  experiments  was 

the  appearence  of  excitations  of  modes  which  were  not 

family  members  of  the  driven  mode. For  example , assume  a  rigid 

cavity  of  dimensions  L  ,L  ,L   is  driven  acoustically  at 
J  x '  y '  z  y 

frequency  w,the  resultant  pressure  standing  wave  is  of  the 

form 

cosk  x  cosk  y  cosk  z  cos(wt+8)  (2.1) 

x       yJ  z 

where   k  =NTT/L   ,  k  =M1T/L   ,  k  =LTT/L  (2.2) 

x      x  '   y      y  '   z      z 

and  N,M,L  are  integers . Eq. ( 2 . 1)  can  be  represented  by  the 
notational  shorthand 

(M,M,L/w,9) 
If  the  cavity  is  driven  to  excite  the  (0,M,0)  mode, then  the 
family  of  standing  waves  consist  of  all  of  those  of  the  form 
( Q ,nM,0 /nw, 9  ) when  nw=nw~  M  Q. 

As  it  is  stated  in  £3j,TTThe  standing  waves  which  can  be 
excited  in  any  real  cavity  deviate  from  the  predictions  of 
the  linear  wave  equation  with  ideal  boundary  conditions 
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"for  the  following  reasons : 

(a). The  presence  of  boundary-layer  losses  at  the  cavity 
surfaces  yields  a  dispersive  contribution  to  the  wave  equation 

(b) . Geometrical  irregularities  alter  the  effective 
dimensions  of  the  cavity. 

Both  of  these  mechanism  can  be  treated  as  equivalent  as 
long  as  the  shift  in  frequency  are  so  small  that  the  actual 
resonances  are  close  to  the  theoretical  values  resulting 
from  the  classical  model."  These  are  treated  by  assuming 
the  dimensions  are  exact, and  the  apparent  phase  speed  is 
determined  on  that  basis. 

The  resonance  frequency  for  each  standing  wave  is 

defined  as  [ 3J 

w  =C  [(n   k  )2  +  (n   k  )2  +  (n   k  )211/2  (2.3) 

n   nL   x   x       yy       z   z  J 

where  k's  are  given  by  Eq . ( 2 . 2 )  and  n  is  a  shorthand  for 

the  set  (n  ,n  ,n  )  where  n  ,n  ,n   are  integers, and  C   is 
x '  y '  z         x '  y '  z  °    '      n 

the  apparent  phase  speed  appropriate  for  that  frequency. 

The  non-linear  wave  equation  applicable   to  this 
problem  can  be  obtained  as  follows : 

The  continuity  equation  for  wave  propagation  in  Eulerian 
coordinates  is 

v.(/a)  +  |f=o  (2-IO 

this  equation  can  be  written  in  terms  of  the  condensation, 

V.C(  i+s)uj  +|5-  =o  (2.5) 

The  equation  of  motion  in  Eulerian  coordinates  for  a 
contained  viscous  fluid  is 
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ridt  C2,6L1 


where 


ps4C11(^/=4k2[,  +  rs  +  ii^ii52J 


ODAT=Other  dispersive  and  absorptive  terms  arising  from 

boundary  effects . 

Eq.(2.6)  can  be  rearrange  in  the  form  of 

|£+(u.<7)a+^vP=^--£«-  (2.8) 

where  the  operator  Lx    describes  those  physical  processes 
leading  to  absorption  and  dispersion. 

One  can  write  V*U=0  and  therefore  ULs^Q  , where  (p  is  the 
velocity  potential , based  on  the  irrotational  velocity 
assumption . Hence  ,  £ u.s  L  V$  .Replacing  U  =  V^)  and  using  the 
condensation ,Eq .( 2 . 8 )  can  be  written  as 

Now, with  the  help  of  Eq.(2.5)  and (2. 9),  and  after  a  good 
deal  of  manipulation  the  non-linear  wave  equation  may  be 
approximated   in  terms  of  velocity  potential 

where 

It  should  be  noted  that  if  the  fluid  is  lossless  and 

c  =Cn  then  (2.1Q)  reduces  to  a  previously  known  non-linear 
p   U  . 

wave  equation  C  9  3 
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c:o^£[cv$/+^(#)2] 

(2.10a) 
In  order  to  express  the  approximate  non-linear  wave 
equation  in  terms  of  acoustic  pressure  and  particle  velocity, 
one  can  rearrange  the  Eq.C2.9)  in  terms  of  p  and  u  and 
combine  that  equation  with  (2. 10). The  result  is  a  quadratically 

non-linear  wave  equation  t  2] 

-a 


(2.11) 


y 


If  it  chances  to  be  that  Vj^iJ   and  \.-F~)     nearly  satisf 

2   2. 

the  wave  equation,  C0  D  (     )  =  ^   ,then  on  the  RHS  of  Eq . 

(2.11)  C072=-2L   and  (2.11)  becomes 

tel&)i-&t*?(7&f+(£fl  ^      i   (2.12) 

Further, if  it  happens  that  jH  f  £   \l  —  .11  ( JL  )     , 

as  is  true  for  solutions  to  the  wave  equation  separated  in 
cartesian  coordinates,  then  [2 J 

^^VoCo1/ -    5?  1  -r  (jfcjJ  J  (2.i3) 

(Note  that  this  is  true  only  for  cartesian  coordinates.) 

As  it  is  stated  in  [3] "The  LHS  of  Eq.(2.13)  is  the 
classical , linear  wave  equation   pertinent  to  the  system 
under  study. The  RHS  can  be  interpreted  as  a  forcing  function 
consisting  of  a  three-dimensional  spatial  distribution  of 
phase-cohorent  sources. In  a  second-order  perturbation  theory, 
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"this  forcing  function  is  obtained  from  the  classical  (first- 
order)  solution  of  the  acoustic  problem. The  second-order 
perturbation  solution  describes  the  nonlinearities  resulting 
from  the  self  interaction  of  the  classical  solution. Higher- 
order  perturbation  solutions  consider  the  interaction  of  the 
non-linear  solution  with  itself, and  the  forcing  function  is 
composed  of  products  of  both  classical  and  nonlinearly 
generated  terms . 

Thus, if  a  system  is  driven  at  frequency  w,the  nonlinear 
term  in. .. "equation  (2.13)  "...  will  force  the  existence  of 
all  integer  multiples  nw  of  the  driving  frequency  and  the 
full  solution  must  contain  all  harmonics  of  the  input 
frequency. In  a  closed  cavity, each  of  those  nonlinearly 
generated  waves  whose  frequency  lies  close  to  the  resonance 
frequency  of  a  standing  wave  of  the  cavity  and  whose 
associated  spatial  function  matches  that  of  the  standing 
wave  can  be  strongly  excited. Just  how  strongly  will  depend 
on   the  quality  factor  Q  for  the  particular  resonance  and 
the  difference  between  the  resonance  frequency  of  the 
standing  wave  and  the  harmonic  nw 

"Consider  two  limiting  cases . 

"(1)  If  the  forcing  function  does  not  have  its  frequency 
nw  close  to  w  ,this  standing  wave  is  being  forced  at  a 
frequency  far  removed  from  its  resonance . Since  this  yields 


the  inequality 

|Co2o2f|»  |£.*H 


(2  .14) 

losses  can  be  ignored  in. . . "Eq. (2 . 13 ) . 
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"(2)  if  nwv  w  ,  then  the  standing  wave  is  being  forced 
near  resonance , and  losses  must  be  retained  in . . . "Eq . ( 2 . 13 ) . 
"The  value  of  C   can  be  determined  from  the  apparent  dimensions 
of  the  cavity  and  the  measured  resonance  frequency  w  . 

"The  losses  are  described  by  the  measured  Q   of  the 
resonance . This  means  that  the  linear-wave  equation  operator 
for  the  system  can  be  written  as 

C0°  *4tL    =  CnD  "  Qn  *  (2.15) 

"Comparison;  of  cases  (1)  and  (2)  reveals  that  the 

responce  of  the  cavity  when  nwjOw   is  order  of  1/Q  compared 

to  that  when  nwvw  .Thus  for  the  high-Q  resonances  usually 

n 

encountered  in  cavities  with  rigid  walls, the  components  in 
the  forcing  function  which  excite  standing  waves  far  from 
resonance  can  be  ignored  compared  to  those  components 
exciting  standing  waves  near  resonance.... 

"The  non-linear , coupled , transcendental  equation  appliciable 
to  this  problem  can  be  expressed  as 


.) 

] 


(cos )  p  n-1  f cos  ? 

sin](/n-9n)=N0M0^n    COsenL1/2    ^x      Vn-jjsinJ    (^<- 

-f>n  +  j    Rj[sin]    'WjO  (2.16) 

for  all  n>l.The  values  of  Q   and  w   must  be  determined  from 

•n       n 

the  infinitesimal-amplitude  behavior  of  the  cavity. The  Mach 
number  Mn  and  driving  frequency  w  are  known  and  Nq  has  the  value 
Nn=l/2  for  a  one-dimensional  standing  wave 
1/M-  for  a  two-dimensional  standing  wave 
1/8  for  a  three-dimensional  standing  wave." 
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r   is  the  Fourier  coefficient  of  nth  harmonic  component, 

normalized  such  that  R..=l.  /   is  the  phase  angle  of  the 

nth  harmonic  component , where  /-.=0.,and  the  phase  angle  9 

is  given  byT3] 

tan9  =-F  (2.17) 

n    n 

where    Fn=Q    [(nw)2-w    2]/(nw)2  =  2Q    (HHzSin)    f  or  iiSlSp  «1      (2.13) 


n  xnw       n  ^n   w 


n 


wn 


Equation  (2.16)  can  be  solved  by  a  method  of  successive 
approximations  on  a  digital   computer . This  has  been  done  by 
[3]  and  [5]  and  both  decided  that  the  theoretical  model  can 
be  used  to  identify  the  modes  of  a  non-ideal , rigid-walled 

cavity  provided  quantities  e   to  be  defined  later  are 

J    r  n 

sufficiently  small. The  theoretical  model  in   its  present 
form  fails  to  account  for  the  excitation  of  modes  other 
than  those  belonging  to  the  family  of  the  driven  mode. 
This  excitation  was  observed  to  occur  only  in  the  case  of 
nearly  degenerate  modes. It  is  believed  to  be  caused  by 
some  linear  coupling  mechanism  within  the  cavity. 

The  purpose  of  this  research  is  to  see  if  the  presence 
of  wall  irregularities  can  explain  how  non-family  members 
may  be  strongly  excited, and  to  present  an  example  to  suppor" 
this  theory. 
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3. DEFINITIONS   AND   NOTATIONS 

A. FREQUENCY  PARAMETER 

The  frequency  parameter  is  a  quantity  which  indicates 
the  position  of  the  driving  frequency  relative  to  the 
resonance  frequency , f. ,of  driven  mode  in  terms  of  the  Q, 
of  the  driven  mode. The  frequency  parameter  is  defined  by 

F1=2Q1(f-f1)/f1  (3.1) 

B. STRENGTH  PARAMETER 

The  investigation  of  the  pressure  waveform  in  the  cavity 
required  the  calculation  of  the  strength  parameter  from  the 
observable  quantities .The  strength  parameter  is  defined  as 

STRPM  =  MQfQ1  (3.2) 

where 

M   -p    / ( p  r    2\  (33) 

u0   1   '0  0  ;  u*  J 

and  P,  is  the  peak  amplitude  of  p, ,the  pressure  distribution 
of  the  driven  mode . 

In  terms  of  observable  or  calculable  quantities  it  is 
reformulated  asT5j 

STRPM=vf2'  V^Q1/(Sm^CQ)2  (3.4) 

where  V  and  S   are  the  rms  voltage  reading  and  microphone 
sensitivity  respectively  of  the  receiver  used  to  sense  the 
standing  wave. 
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C.  e   is  defined  to  indicate  the  position  of  w   relative 
n  r  n 


to  the  classical  harmonic  frequencies ,nw, . 


w  -nw-i 
e  =  J2 i 

n     nwi 


(3.5) 


and  one  can  relate  e   with  F   from  (2.18)  such  as 

n       n 

F  =  2Q  e 
n   xn  n 


C3.6) 


D.  A  pictorial  representation  of  F  which  will  be  useful 
throughout  the  development  is  given  in  Fig. 1. From  now  on 
three  subscripts  will  be  used  for  convenience,  i.e. 
F   becomes  F 


n 


nml ' 


1_ 
Qnml 


FIGURE   1 


,=[i.c-Ss^)2] 


nml 


nw 


)  JQ   t  =cotf 
J  xnml    wni 


nml 


(3.7) 


0 


nml 


=  S 


Qnml  Sin(^rml""  V.  _2 """"'"  1/2'  Jnml 

nml 


(3.8) 
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4.  THEORETICAL   DEVELOPMENT 


A.  CAVITY  DESCRIPTION 

Assume  a  perfectly  rigid-walled  rectangular  cavity  which 
has  one  wall  perturbed  such  that  the  cavity  dimensions  are 
L  L1+  £ f (y , z)J  ,L   and  L   as  shown  in  Fig. 2  below. Also  assume 
the  perturbation  on  the  boundary  is  very  small  compared  to 
the  cavity  dimensions,  (6f(y,z)l«  1. 


x  Perturbed  wdll 


->^ 


FIGURE  2 
The  cavity  is  to  be  excited  by  a  source  near  the  origin 
in  such  a  way  that  the  (M,M,L)  mode  is  driven  at  a  frequency 
close  to  its  resonance  frequency. 
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B.  THE  PERTURBED  BOUNDARY 

For  a  rigid-walled  rectangular  cavity  with  ideal  boundaries 
(€=0),the  pressure  p~  obtained  from  the  linear  wave  equation 
with  losses 

is  subject  to  the  following  conditions, 

Vpn.n=0      at  x=0,L  (4.2) 

y  =  0  ,L 

J        y 

z  =  0  ,L 
'  z 

where  n  is  the  local  normal  to  the  ideal  boundary. The  solution 
for  pn  in  terms  of  Mach  number  is  given  by  [4] 

p0 
„=Mn  cosk  x  cosk  y  cosk  z  cos(wt+9) 

/°0c02    °       x        y  (4.3) 

=MQ  (N,M,L/w,9) 

where  k's  are  given  in  Eq . ( 2 . 2 ) ,  and 

w  =  c   (k  2+k  2+k  2)1/2  (14.4) 

p    x    y    z 

If  the  cavity  has  perturbed  walls, the  solution  will  be 
in  terms  of  a  summation  of  the  classical  linear  solution 
for  ideal  boundaries  plus  perturbation  correction  terms 
due  to  the  irregular  boundary: 

p=p0+€p'+62p' '+ (4.5) 

Since  the  magnitude  of  the  boundary  perturbation  is  kept 
small, second  and  higher  terms  in  £  can  be  considered 
insignificant , so  that 

p=pn+€p*    (to  first  order)  (4.6) 

and  p  must  satisfy  the  following  conditions, 

□  2p  =  0  (4.7) 
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and 


Vp.n=0 


at  x=0,L  [l+€f(y,z)] 


(4.8) 


y=0,L 

z  =  0  ,L 


y 


where  n,the  local  normal  to  the  real  surf ace, is  obtained  by 
taking  the  gradient  of  the  equation  for  the  boundary , given 
by  [2] 


n 


=  v[x-Lx[n-€f(y,z)]} 


Thus, to  the  first  order  in  £, 

n-X-fcL  — r y-feL   — pr —  z 

and  when  Eg. (4.11)  is  used  in  (4.8)  the  result  is 


(4.10) 


(4.11) 


f3p    cL  ?SiL>i2:.eu'5f^22  2tl  -0 


(4.12) 
x[»+£fC**)] 

A   Taylor    series    expansion [4]  for   p    evaluated    at    the    real 


boundary    L    [l+Gf  (y  ,  z  )]    produces 

ol  -  P I  +  — 


£Lxf  C*,2> 


x  =  L: 


+  }M       CeUf^D]- 


2  0x~    x,Lx 


(4.13) 


Substituting  Eq . ( 4 . 6 )  into  RHS  of  ( 4  .  13 ) , taking  the  partial 
derivative  with  respect  to  x  on  both  sides  and  keeping  the 


first-order  terms  in  €  , yields 
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Ix. 


+  £ 


Df' 


Ixsl* 


X=U  ^  'x*Lx    (4.14) 


'a*  L_,_  1& 


Taking  the  partial  derivatives  with  respect  to  y  and  z  and 
using  exactly  the  same  procedure  gives 


If. 

^    'x«Lxti+cf<y.D3 


+  £ 


U-Ly 


2fl 


X--L; 


eLxf^2\r:i'5) 

X"-L* 
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,2£ 

x--L*0+£f<*"J9z 


+  G 


3f' 


X^L*   2)2 


Jr-lv 


(4.15) 
Substituting  ( 4 . 14 ) , (4 . 15 )  and  (4.16)  into  (4.12)  and  keeping 
the  first-order  terms  in  £.  results  in 


It 
3X 


,»l/  *L~   '   W    ^    ^    ^   ^2  JX=LX  ( 


X=LV  *-"•*  (4.17) 

The  RHS  of  Eq.(4.17)  can  be  represented  as  a  Fourier 
series  in  cosines, so  that  p'  can  be  expressed  as  a  summation 
of  normal  modes. Hence, 

If.' I        =  J  f  amlcos^ijcosHz  cos(wt+9) 

0A~    'XiLx       m.Q    l=0 


or 


(4.18) 


2*1 


oo       co 


=  ££dmt(0,m,l/w,9; 

X-L*.         nr\i0     UO 


where 


*«  -±1/"'* 


dz 


A.  -JtJJ  [««¥»* 


2.     r 


0    0 

Ly  Li 


Lot      - 


L^L2  , 


j-G]  Co*i£   djdz 


0    0 


(4.19a) 


(4.19b) 


(4.19c) 


(4.19d) 


and   G    is    defined    as 


6a¥ 


x=L 


(4.19e) 
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In  order  to  find  the  contribution  to  p'  from  each  of  the 
terms, it  is  stated[2] that  for  a  cavity  forced  by  a  dynamic 
boundary  condition  at  a  boundary 


2   2  . 

Ql  p 


=  0 


(4.20) 


and  the  dynamic  boundary  condition  from  the  (m,l)th  term  is 


3x 

then 
i 


=  Acos^j  cosiUz.   e 


i(wt+0) 


X*L, 


ao 


*    -A(tf^^ 


X>«c-^«E*"*i;*"srxie 


(4.21) 


iW+Q+<Jnm^ 


where 


=f 


n  /  2 


if  n  =  0 

if  n=l,2,3, . . . 


(4.22) 


(4.23) 


and  S   -  is  given  by  Eq.(3.8). 

Applying  this  solution  to  Eq.(4.18)  gives  the  complete 

solution  for  the  first-order  perturbation 
oo 

f  =  ZL/Pml 

00  (4.24) 


rruo 


----±— r-  &LM-i)nSnml  (n.m.L/w^  +  Cin) 


and  combining  (4.24)  with  (4.6)  yields  the   total  acoustic 
pressure  in  the  cavity. 


f  =  fN,M,L/w,6)--^r-  £fWM-u  W"''"'1-  Me+<w) 

1=0 

If  this  is  near  a  resonance ,w^v  w   .  ,then  this  term  will 

nml 

dominate  the  summation  and  all  other  non-degenerate  terms 
can  be  omitted. Consequently , the  Eq.(4.25)  becomes 

f«(HM.L/^0— S-T  ^LAa(-l)nSnmL(n,m(L/^e^07mi)(4  26) 
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5.  SPECIFIC  EXAMPLES 


A.  CAVITY  WITH  STEPPED  PERTURBATION 

Assume  that  the  rigid-walled  rectangular  cavity  given  in  Fig. 2 
is  perturbed  as  shown  in  Fig. 3  below, and  also  assume  that  the  cavity 
is  driven  in  the  (0,1,0)  mode  resulting  (0,l,0/w,9)  standing  wave 
and  that  the  (0,2,0)  and  (1,0,0)  modes  are  (nearly)  degenerate. 


•mic. 

-> 

A 

<- 
L 

\L 

.source 

FIGURE  3 
From  Fig. 3  the  equation  for  the  boundary  at  L^  can  be  found, 

By  means  of  Eq.(4.8)  £  and  f(v,z)  can  be  written  as 

f(j,Z)--Ci(!*-L')  -  KM-)]      it     lokl 

Since  the  (0,2,0)  and  (1,0,0)  modes  are  degenerate  the  emphasis 


(5A.1) 

(5A.2) 
(5A.3) 


of  this  development  will  be  on  these  particular  modes. The  pressure 
distribution  of  (0,2,0)  mode  is 

?020  =  P2  cos^  cos  (2*i  +  9z) 

or  (5A.4) 

-Po2o=P2  (O,2;O/2W,02) 

where  P„  is  the  amplitude  of  (0,2,0)  mode. 


Utilizing  the  theory  developed  in  the  proceeding  sections  and 

using  the  equations  (4.17)  through  (4.26),  the  first-order  perturbation 

correction  can  be  found 

9^,z)   B'Poao 


-uc 


XsL, 


^ 


oy 


L.L 


X*  Lie 


(5A.5) 


ax 


Eq.(5A.5)  can  be  written  as  a  Fourier  series 


(5A.6) 


m=0 


Inversion  of  the  Eq. (5A-5)  and  (5A.6)  yields  the  Fourier  coefficients 

3m=^,[5io(flU')c0S(^L')-Si',(^gC0S(Jffi!L)J  (5A.7) 

and 

30  ,d.[s,n(^L')-^(^L)]  (5A8) 

Recalling  Eq.O+,23)  and  (4.24) , first  order  perturbation  correction ,p' , 


is  found  as 


LaJdE  j£_  (2)(-i)a,oo5,n(3roo  c°*x 


jC2wt+e2+*r0o)         g) 


and 

€f  =  ip6a0JLria,00*«"0ii.cosJg.  e  (SA.io) 


where 


(gN^f-C^f 


Hence, the  total  pressure, associated  with  the  angular  frequency  2w, 
in  the  cavity  is 

^  Pz(o,2, o/z#,Bz)+  ^k.Ga0q,0O5in<^0  (i,o,o/2u/,et+<s0) 

(5A.11) 


25 


The  total  pressure  at  the  microphone  oosition,x=0  and  y=L  ,is 

i(2wt+e2)     .  i(2u/-t+e2  +  c?w) 

(5A.12) 
'mic.  position 


1  'mic.posi'+ton  s"  J 


Define  B=[        J  and  after  a  little  manipulation  and  use  of  trigonometric 
identities , ( 5A. 12 )  becomes 

f  LicposWdT  P2|C«  +  BCoS<^0)Cos(2wUea)-(Bsin<?Too)  sm(2wt-,92)j 

(5A.13) 
and  the  amplitude  of  the  total  pressure  in  the  cavity  is 

P  L.p<*iHon  ■  P2\/0  +  B<:°Scroo)2+(e>S1nC77M)1  (SA.UO 

Eq.(5A.lM-)   is  the  corrected  value  o^  the  amDlitude  of  the  second 
harmonic  of  the  driving  mode, obtained  bv  Eq.  (2.16) , because  of  the 
boundarv  irregularity  given  in  Fig. 3. 

New, it  is  desired  to  express     sin<3f]J0  in  terms  of  the  frequency 
parameter  of  the  driving  mode, (0,1,0)  .With  the  help  of     Fig.l,sin0f^o 
can  be  written  as 


'  HOC    ' 

2f 


{i*«ft.ci-(Wrj: 


If  f  approaches  to  zero  then  <7^0  approaches  to  "IT  ,and  if  f 
approaches  to  infinity  then  (7^j0  is  close  to  zero. In  these  same 
limits  cos  Cj~  goes  to  -1  and  +1  respectively.  Hence 

COS  C700  ^ijl-sin1^  (5A.16) 
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For     2f  "sj  f1QQ      ,(5A.15)  becomes 

Sm^oo  *  c r; JF-TTzT'h  C5A.17) 

*"  TiOO  J 

Recalling;  Eq. (3.1)  and  (3.5),F1Q0  and  e10Q  can  be  written  in  the  form  of 

FblO   s  2C^0.O  ^-^  (5A.18) 

and 

p     _  f\0o  -  2ft>to  (5A.19) 

Eq.(5A.19)  can  be  solved  -For  finn  and  this  substituted  into  (5A.17) 

^^too  -7 ^010  —    -_ _- _  (5A.20) 

rioo  Woio       f0,0  [\  +  eioo) 

Use  of  l/(l+einn)=l-einn  and  little  manipulation  reveals 

awOSc   =        ■  i  (5A.21) 

v1+[|^^F--2^'»e»»o]Ci-e,o,)]z 

As  a  result, equations  (5A.1U) ,  (5A.16) ,  (5A.21)  are  the  final  amplitude 
correction  of  the  second  harmonic  o^  the  driving;  mode  obtained  bv  Eq .  (2 .  16 ) 
A  computer  program  for  this  was  develooed  bv  author  and  is  given  in 
appendix  A. 
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B.   LINEARLY  PERTURBED  CAVITY 

Using  the  same  assumptions  in  section  A, assume  that  the  rigid-walled 
rectangular  cavity  is  perturbed  linearly  as  shown  in  Fig. 4-  below. 

►Wc  y 


ih* 


Saurtt 


\ 


\ 

U-JL. 


FIGURE  4 
The  equation  for  this  perturbation  is 

x.L„  [■+£.(!-.£»)] 

Hence , 


and 


6-f 


(5B.1) 


(5B.2) 


(5B.3) 

Applying  the  same  procedure  as  in  section  A, the  Pirst-order  perturbation 
correction  and  the  total  acoustic  pressure  associated  with  angular  frequency 
2w  in  that  oarticular  cavity  can  be  found 


3x 


.'I     .  Lj[-J^**  «*  cos(2Wtt^)](-f  )] 


=  !gAb  sirt-*^  Cos(2wl+  92) 


(5B.4) 


12, 

Eq.(5B.<+)  can  be  written  as  a  Fourier  series  and  the  Fourier  coefficient, 
a  ,is  found  bv  an  integration  Drocedure  evaluated  in  the  interval 
0  to  L  .The  result  is 

y 

a  s.l{   C0S(2-m)TT^   C03a4m)Tr J L.1,^ 

^  IT  I       (2-m)  C2-fm)  (2-m)     (2-fm)  / 


c9o      =  0.0 
3z     =  0-0 


(5B.5) 
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Recalling  the  Eq, CiK 2 41, the  first-order  perturbation 
correction  for  CQ.,  2,0.1  mode  is 

f '  =  _  iE^u  ^_  c2)c-j)1510o  (i.t>.o/2v/,el*a-u>0 )      (5B.6) 

and  the  total  acoustic  pressure  associated  with  angular 
frequency  2w  in  the  cavity  becomes 

f  =  (0<2,0/2uJ,92)-^4^2r2)("1)i3,oo^'0'0/2W'924^0o)  (5B'7) 
since  a~=0  .  0 

According  to  the  calculation  developed  above  there  is 
no  need  to  make  a  first-order  perturbation  correction  to  the 
CO, 2,0)  mode  in  the  cavity  shown  in  Fig, 4,  As  a  result, the 
pressure  distribution  is  equal  to  the  second  harmonic  of  the 
driven  mode, since  an=Q.O  and  this  yields  «p'=Q.Q 
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C.  CAVITY   WITH   WEDGED   PERTURBATION 
Under  the  same  assumptions  made  in  sections  A  and  B, 
assume  that  the  cavity  is  perturbed  as  shown  in  Fig. 5  below. 


ih 


mic 


»5oi)ree 


FIGURE  5 


The  equation  for  this  perturbation  is 
By  means  of   Eq.(4.8)  £  and  fCy,z)  can  be  written 


(5C.1) 


as 


(5C.2) 


and 


(5C.3) 


f(*z)  =0.0  >  X  4.   L^/2 

Applying  exactly  the  same  procedure  followed  in  section  A, 
the  total  pressure  amplitude  in  the  cavity (in  first-order 
perturbation)  is 


=   p2VCl+Bcos(^o)Z+C65in(7700) 


mic.  posifibo 


(5C.4) 


where 


5    -  -L-<3o€Q,ooSJnar00 
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sinCJ^0  and  cosiyfJ0  are  given  by  Eq.(5A.21)  and  C5A.16)  respectively. 

The  theoretical  predictions  of  these  specific  examples 
were  examined  with  series  of  experiments  developed  by  [8  J  . 
The  further  discussions  about  these  will  be  given  in  the  next 
section. 
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6.  RESULTS 

In  this  section  the  theoretical  predictions  performed 
in  sections  5A,  5B  and  5C  will  be  compared  to  the  experimental 
results  obtained  by  £8l . 

The  information  on  the  empirical  losses  and  resonance 
frequencies  is  contained  in  the  Q's  and  e's. These  are  the 
values  used  in  the  computer  program  to  predict  the  harmonic 
distortion  on  the  basis  of  Eq.(2.16)  and  is  plotted  as  thin 
solid  curves. The  results  of  including  the  first-order 
perturbation  correction  are  plotted  as  thick  solid  curves 
for  each  specific  example. The  theoretical  curves  in  figures 
6  through  16  were  plotted  along  with  the  experimentally- 
measured  values  for  the  cavity  configurations  shown  on  top 
of  each  figure. The  theoretically-predicted  values  were 
generated  for  frequency-parameter  intervals  of  0.2,  and 
the  experimentally-measured  values  were  plotted  as  square- 
blocks. Data  were  taken,  and  theoretical  predictions  made, 
for  different  strength  parameters  for  the  (0,n,0)  mode 
associated  with  different  cavity  configurations.lt  is 
important  to  note  that  the  n=2  distortion  peaks  when  the 
system  is  driven  at  this  frequency .That  is,  when  the  driving 
frequency  w  is  equal  to  (1/2)  w~ , there  is  maximum  content 

of  P-.The  point  where  this  occurs  for  each  P  /P-  curve  is 
2  n   1 

indicated  byyyarrow  with  the  label  of  F  _n.At  that  point 
the  value  of  frequency  parameter  is 

F020=2  Q010  e020 
(The  same  thing  could  be  done,  of  course,  for  any  member 
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of  the  (0,n,0)  family).  The  arrow  labelled  as  Finf  indicates 
the  position  where  the  nearly  degenarate  (1,0,0)  mode  is 
resonant,  and  the  value  of  Finn  is 

Fioo=2  Q010  e100 

The  theoretical  predictions  made  in  section  5B  were 

compared  to  the  experimental  results  as  seen  in  Fig. 7. 

When  Fig. 7  is  compared  to  Fig. 6,  the  unperturbed  cavity, 

it  is  clearly  seen  that  the  theory  and  experiment  are 

excellently  in  agreement. 

For  a  wedged  perturbation,  the  theory  predicts  the 
frequency  of  the  second  harmonic  at  which  the  effect  of  the 
perturbation  occurs  as  seen  in  Fig. 8. The  predicted  magnitude 
of  the  perturbation  effect  for  this  configuration  is  in  good 
agreement  with  the  experiment .The  anomolous  behavior  of  the 
third  harmonic  in  Fig. 3  is  unexplained. 

For  stepped  perturbation,  when  the  cavity  is  perturbed, 
but  the  geometry  leads  to  no  predicted  correction  as  seen  in 
Fig. 9  or  leads  to  predicted  correction  less  than  about  0.02 
as  in  Fig. 12  or  less  than  about  0.05  as  in  Fig. 14,  then  it 
was  observed  that  there  was  very  little  or  no  effect  from 
the  (1,0,0)  mode .Agreement  for  these  cases  is  good  except 
for  the  region  lying  between  frequency  parameter  4  and  9 
in  Fig. 9. What  happened  in  that  region  is  also  unexplained, 
but  it  was  observed  one  time  only. When  the  amount  of 
perturbation  correction  is  increased  the  theory  predicts 
effects  larger  than  experimentally  observed .However ,  the 
effect  of  the  perturbation  appears  at  the  right  frequency 
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parameter  as  is  seen  in  Fig.'s  10,11,13,15  and  16. 

Choosing  the  shim  position,  length  and  magnitude  is  very 
important  as  well  as  is  choice  of  the  strength  parameter. 
For  the  shims,  A  =0.04-  and  0.25  inches  for  stepped  and  wedged 
perturbations  respectively .The  effect  of  strength  parameter 
can  be  seen  in  Fig.'s  15  and  16. The  experimental  data  associated 
with  the  third  harmonic  in  Fig. 15  were  believed  to  come  from 
harmonic  distortion  in  the  piston  motion. 


33 


1,0  [MOC£  :QJQ.   rQ(0":2?fl.^6TE(nrQ.Q       ]Q|QQ:3S7" 
5TfiP,vi:p.59Q!Qm:5.qi.22  'EC2):7.9g,D^r 


0.00 1 


-2-1  0  1  2 

FREQUENCY        PARAMETER 


FIGURE    6 


34 


1.0 


MODE'.QIO     |Qfn=  ^5316  rECO=Q-0      !QiOO  =304-35 


STRPM:0.365iQc2)-.^o2.5   £C2)=80.7«IO 
NTQP:2d    ;Q(.5)»t/VB.9TE(7)=8S.5«ro5~' 


r-fc-8- 


.(4)^05.5  ;E(4)=15.^'CV| 


E 100;  2.22  xm 


at 
skL. 


-THSPftV  1 


□  EXPERIMENT 


FREQUENCY  PARAMETER 

FIGURE    7 


35 


1.0 


MODE  L-QJLQKU j )  =2  S3J7rQJoo£[£^K 
STRPM  :Q222  T CU  2\=  550.7! 


0.001 


E-  ioorlo.i<,io 


E(0=  o-Q 
E(Z)s-7.S»uf 


IOEAL. 


^1  .zlj— '  PEgTilHaEp 
EEXPE21MENT 


FREQUENCY        PARAMETER 


FIGURE    8 


36 


CL-IQ-- 


37 


CL   IQ_- 


38 


"s'iPM  : 0305  QfZ3  '■  377         "  EC3 : ZJ «.  5 1 


nTcP;20       I Q  (5V1*  SO 
TsTMAx  :  iO- 


JD6At_ ] 


£^5)'.  t.+  xic"-  —      '         G  EXPEeiJr^T 


:.:  : 


FREdUENa  F  .  -     v  ETER 


39 


1.0 


MODE '.  OIO    \Qi 0 :.  iSfc.57 Gj  IQO ^3^751  £ ]0O^±f}'^  EH) :_0. 0 


E(2)-.-6,^idJ;  ^ 


T  I —  IDEAL 

5.9*2  —  PE£TUC3EC 


0.00  Ik 


2  3 

FREQUENCY       PARAMETER 


FIGURE    12 


40 


FREQUENCY        PARAMETER 
FIGURE    1 3 


41 


-3- 


Moce :  oio    { q < i  v  2 fe.7.^  JQ  <ooT572rj~e>op  ^~8*~fc?~E( *> :_ a o 


T IDEAL        | 

—  P6«TUftBClj> 


0.1 


OJOI 


OCOI 


FREQUENCY         PARAMETER 
FIGURE    14 


42 


1.0 


MQOEi0iCL_jQ(«k2  78-bj.  'Qioo-369.96 

SJRPMr  0,205  ;  Q12L-.38ZJL  J 

NTQP.2Q      |Q(^)^68.53 

NMAX  -.   IO 


—  ^>VTC 


£  IQO'-TfrxIcTjE  ( I )  =  jD^O 

~£i2)»  20,6*10 


O 

HI  <(V 


f  j IDEAL 

6"!—  PERTURBED 


□  EXPERIMENT 


FREQUENCY  PARAMETER 

FIGURE    15 


'•°  [mode  :  oio]  q(0  _:276.~6S  Q  100  t3fe7.gTeTb°^-8»'Q  E'O    °°     -T 
^TiiPM  :  0.29    Q  (2)  :  580;r  E(2^  :  50.2«io 

NTop:20      |q(3):V7\5.7 

tvlMAX  :  !0 


01 


OJO  i 


0.001 


-^.i-jir 

;*■ 

■  < 

■  •          '        • 



...U44...;...;...!_.U..| 

...  :_ 

£3£ESJ 

EfSE^i 

^«J 

.Vtt 


/ 


15 


7* 


FrZ0i 


FREaUENCY         PARAMETER 


0" 


4 


FIGURE    16 


44 


7.  CONCLUSION 

Non-linear  theory  has  been  applied  to  standing  waves  in  a 
rigid-walled  rectangular  cavity  with  a  perturbed  boundary 
in  order  to  find  one  possible  mechanism  for  the  excitation 
of  a  standing  wave  other  than  those  belonging  the  family  of  the 
driven  mode . It  was  observed  that  such  an  excitation  exists 
if  the  boundary  perturbation  and  the  dimensions  of  the  cavity 
are  favorably  chosen. 

It  appears  that  the  present  theoretical  model  succesfully 
predicts  the  major  features  of  harmonic  content  for  finite- 
amplitude  standing  waves  in  the  cavity  when  the  geometry  leads 
to  no  perturbation  correction  (Fig. 7  and  9). When  the  magnitude 
of  the  perturbation  is  increased  the  predicted  features  were 
larger  than  experimentally  observed . Second-order  perturbation 
corrections  may  be  needed  to  account  for  these  discrepancies. 
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APPENDIX  A 

The  original  computer  program  for  Eq.(2.16)  was  prepared 
by  Coppens  in  197  3,  and  author  made  some  extentions  to  that 
program  so  that  it  would  (1)  calculate  the  perturbation  cor- 
rection and  (2)  present  the  results  graphically.  The  program 
calculates  the  relative  amplitudes  and  phase  angles  of  stan- 
ding waves  in  cavities  keeping  the  strength  parameter  constant 
and  changing  the  frequency  parameter  to  generate  response 
curves  showing  the  amplitudes  of  the  nonlinearly  excited  stan- 
ding waves  as  function  of  the  frequency  parameter  Ej  •  It  also 
calculates  the  perturbation  correction  according  to  Eq.(4.2^) 
and  then  finds  the  total  relative  pressure  amplitude  using 
Eq.C4.26).  The  program  also  draws  the  graph  of  the  relative 
pressure  amplitudes  of  the  ideal  cavity, total  relative  pressure 
amplitude  of  the  perturbed  cavity  and  the  experimental  data 
on  a  3  cycle  semilog  paper. The  Versatec  Graphics  Plotting 
Manual  [IQ]  was  used  for  the  graphical  processes  on  the  IBM 
36 Q  of  the  Randolph  Church  Computer  Center, Naval  Postgraduate 
School . 
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SOME  USEFUL  INFORMATION  ABOUT  COMPUTER  PROGRAM 
*    Quantities  marked  with  (*)  must  be  controlled  or  changed 

for  each  run. 

*K0N        The  number  of  iterations  throughout  the  region 
of  interest. For  this  program  the  iterations  are 
performed  with  0.2  intervals. 

"NCUR  The  number  of  the  experimental  curves  to  be 
drawn  +  1 

"NDAT       The  number  of  experimental  data  in  the  region 
of  interest 

BUR(I,J)  The  array  that  stores  the  experimental  data 

"XL  The  length  of  the  cavity  in  the  x-direction 

*YL  The  length  of  the  cavity  in  the  y-direction 

"DELTA  The  magnitude  of  the  perturbation 

-STRPM  Strength  parameter 

"FREQ        Frequency  parameter  stored  in  ATA(I,1)  and 
ZER( I , 1) . Input  as  the  maximum  value  of  FREQ 
in  the  region  of  interest 

XDAT(JET)    The  x  coordinate 

YDAT(JET)    Value  of  the  curve  f(x)  for  XDAT 

ATA(I,N) ,N*lThe  array  that  stores  the  logarithmic  value 
of  the  pressure  amplitudes  of  the  harmonics 
of  the  driving  mode 

ZER(I ,N) ,N^lThe  array  that  stores  the  linear  value  of  the 
pressure  amplitudes  of  the  harmonics  of  the 
driving  mode 


:Q(I) 


-E(I) 


Quality  factors  of  driving  mode  and  harmonics 
of  it 

e's  value  of  driving  mode  and  harmonics  of  it 
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*Q100         Quality  factor  of  (1,0,0)  mode 

*E100         e-value  of  (1,0,0)  mode 

*XDAT( JET+1)  Integer  value  of  left-hand  corner  on  the  x-axis . 
It  must  have  the  same  value  as  the  7th  argument 
of  subroutine  CALL  AXIS  for  x-axis. 

*YDAT( JET+1)  Integer  value  of  left-hand  corner  on  the  y-axis . 
It  must  have  the  same  value  as  the  7th  argument 
of  subroutine  CALL  AXIS  for  y-axis . 

XDAKJET+2)] 

'  Increment  value  of  x  and  y  for  scalling  purposes 

YDAT( JET+2), 

HUM  The  linear  value  of  the  total  acoustic  pressure 

amplitude  associated  with  angular  frequency  2w. 

*B  -S-Q.nn  sinGTnn ( 1/4 )  L   an  for  stepped  perturbation 

TT  x100      100        y   0         ri    r 

-£-  Q-jqq  sin(JT"nn  (an/2 )  for  wedged  perturbation 
and  an  is  the  0th  Fourier  coefficient. 
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